AEDC-TR -68-224 

(?  ,  U 

M  I 


— — 


f  " 


/? 


Si  1  Ob’J 

may  ?^1969 

0Cf  1  !970 


FREE-FLIGHT  RANGE  TESTS  OF  THE 
APOLLO  COMMAND  MODULE 


W.  R.  Lawrence  and  W.  S.  Norman 
ARO,  Inc. 


January  1969 


This  document  has  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited. 


von  kArmAn  gas  dynamics  facility 

ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 
AIR  FORCE  SYSTEMS  COMMAND 
ARNOLD  AIR  FORCE  STATION ,  TENNESSEE 


PROPERTY  OF  U,  S,  AIR  FORCE 

AFT'C  library 

F4QbOOf.6y.i-G  -.0001 


? 


mim 


When  U.  S.  Government  drawings  specifications,  or  other  data  are  used  for  any  purpose  other  than  a 
definitely  related  Government  procurement  operation,  the  Government  thereby  incurs  no  responsibility 
nor  any  obligation  whatsoever,  and  the  fact  that  the  Government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by  implication 
or  otherwise,  or  in  any  manner  licensing  the  holder  or  any  other  person  or  corporation,  or  conveying 
any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be 
related  thereto. 


Qualified  users  may  obtain  copies  of  this  report  from  the  Defense  Documentation  Center. 


References  to  named  commercial  products  in  this  report  are  not  to  be  considered  in  any  sense  as  an 
endorsement  of  the  product  by  the  United  States  Air  Force  or  the  Government, 


AEDC-TR -68-224 


FREE -FLIGHT  RANGE  TESTS  OF  THE 
APOLLO  COMMAND  MODULE 


W.  R.  Lawrence  and  W.  S.  Norman 
ARO,  Inc. 


This  document  has  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited. 


AEDC-TR-68-224 


FOREWORD 


The  work  reported  herein  was  done  at  the  request  of  the  Arnold 
Engineering  Development  Center  (AEDC),  Air  Force  Systems  Command 
(AFSC),  under  System  920E  in  cooperation  with  the  NASA  Manned  Space¬ 
craft  Center,  Houston,  Texas. 

The  test  results  presented  herein  were  obtained  by  ARO,  Inc.  (a  sub¬ 
sidiary  of  Sverdrup  &  Parcel  and  Associates,  Inc.),  contract  operator  of 
the  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under  Contract 
F40600-69-C-0001.  The  tests  were  conducted  during  the  period  from 
January  25  to  March  10,  1967,  under  ARO  Project  No.  VT0744,  and  the 
manuscript  was  submitted  for  publication  on  September  13,  1968. 

This  technical  report  has  been  reviewed  and  is  approved. 

Eugene  C.  Fletcher 
Lt  Colonel,  USAF 
AF  Representative,  VKF 
Directorate  of  Test 


Roy  R.  Croy,  Jr. 
Colonel,  USAF 
Director  of  Test 
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ABSTRACT 


Results  are  presented  for  free -flight  range  tests  of  the  basic  Apollo 
Command  Module  configuration  (smooth  heat  shield)  conducted  in  Range  G 
at  nominal  Mach  numbers  of  6.  0  and  8.5.  The  Reynolds  number  based  on 
model  diameter  and  free-stream  conditions  varied  from  50,  000  to 
240,  000.  Both  drag  and  stability  data  are  presented  in  addition  to  trim 
angles  of  attack  corresponding  to  three  lateral  center-of-gravity  positions. 
Comparison  of  the  results  with  data  from  wind  tunnels  shows  at  most 
only  small  sting  effects  on  the  wind  tunnel  data. 
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NOMENCLATURE 


A 

Ai,  A 2,  A3 
CL  C2 

CA 

CD 

Cl 

CLq* 

Cm 

Cmu 

Cmpa.  s„ 

Cmq>  Cm# 
CN 
Cn 
cn/3 

cy 

D 

d 


Reference  area  of  model,  7rd^/4 
Constants  (Eq.  (Ill- 2)) 

Constants  in  representation  of  pitching  moment 
(Eqs.  (IV-23)  and  (IV-24)) 

Axial-force  coefficient 

Drag  coefficient 

Lift  coefficient 

Lift-force  derivative 

Pitching -moment  coefficient 

Pitching-moment  derivative 

Magnus  parameters 

Damping-in-pitch  derivatives 
Normal-force  coefficient 
Yawing-moment  coefficient 
Yawing -moment  derivative 
Side-force  coefficient 
Damping  parameter 

Coefficients  of  solution  (see  Eqs.  (IV-52)  and  (IV-53)) 
Diameter  of  model 

Amount  of  center -of-gravity  offset  (Fig.  2) 
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F,  G 

Fl>  f2 
Ix>  Iy»  Iz 


i 


i,  h  k 


ii>  h>  ki 


Inertia  parameters,  defined  by  Eqs.  (IV- 3)  and  (IV-4) 
Defined  by  Eq.  (IV-44) 

Moments  of  inertia 

Imaginary  unit,  \  -1 

Unit  vectors  in  body-fixed  system 

Unit  vectors  in  nonrolling  body  axis  system 


ir,  jr,  kj.  Unit  vectors  in  range  system 

Ki,  K2,  K3  Defined  by  Eqs.  (IV-33),  (IV-34),  and  (IV-38) 


ka 

kt 

M 

M0 

m 

N 

P,  Q,  R 


q 

Reood 


Re2d 

S 

Tl,  T2 
t 

U,  V,  w 
Ui,  Vi,  wx 

V. 

X 

x 


V  Ix/md2 

\[  Iy/md2 


Pitching  moment 
Mach  number 
Mass  of  model 
Yawing  moment 

Roll  rate,  pitch  rate,  and  yaw  rate  (usually  radians /ft) 


(e/d) 

0. 03384 


(see  Appendix  IV) 


Dynamic  pressure 

Reynolds  number  based  on  free-stream  conditions  and 
diameter  of  model 

Reynolds  number  based  on  conditions  behind  normal  shock 
and  diameter  of  model 

Defined  by  Eq.  (IV-51) 

Defined  by  Eqs.  (IV-54)  and  (IV-55) 

Time 

Velocity  components  relative  to  body-fixed  axis  system 

Velocity  components  relative  to  nonrolling  body  axis  system 

Velocity  of  model 

Axial  distance  along  range 

Axial  location  of  center  of  gravity  (Fig.  2) 
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x,  y,  z 

xi,  yi,  i\ 
xr»  Yr>  zr 
Y 

y  +  iz 

Z 

Zl 


«t 

«1 

H 

Hr 

Hi 

y 

I 

ba>  6j3 

e 

0i 

Ml,  M 2 
? 

P 

^  2 
wi*  w2 
“A.  WB 


Body-fixed  axis  system 
Nonrolling  body  axis  system 
Range  axis  system 
Side  force 

Transverse  displacement  of  model  from  range  axis 
(see  Eq.  (Ill - 6)) 

Force  in  z  direction 

Defined  by  Eq.  (IV-50) 

Angle  of  attack  in  body-fixed  axis  system 
Total  angle  of  attack 

Angle  of  attack  in  nonrolling  axis  system 
Trim  angle 

Trim  angle  for  zero  roll  rate 

a  -  at 

Angle  of  sideslip  in  body-fixed  axis  system 
Angle  of  sideslip  in  nonrolling  axis  system 
Yaw  angle  in  nonrolling  body  axis  system 
Cosine  ae 

Mean  total  angle  of  attack  (see  Eq.  ( III  -  9 ) ) 

Phase  shift  angles  (Eqs.  (IV-52)  and  (IV-53)) 

Elevation  angle 

Elevation  angle  in  nonrolling  body  axis  system 
See  Eq.  (HI-2) 

Dependent  variable  in  the  yaw  equation  (see  Appendix  III) 

Range  air  density 

Roll  angle 

Initial  roll  angle 

See  Eq.  (III-2) 

See  Eq.  (IV-49) 

Defined  by  Eq.  (IV-45) 
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w a,  Wj3  Defined  by  Eqs.  (IV- 30)  and  (IV- 37) 

(  )'  =  d(  )/dx 


SUBSCRIPTS 

o,  t  Trim  conditions 

eg'  Center  of  gravity 

Note:  The  coordinate  system  used  in  this  report  has  the  x  axis  in  the 
direction  of  the  velocity  vector,  whereas  the  coordinate  system  used  in 
Ref.  1  has  the  x-axis  opposite  to  the  velocity  vector.  Therefore,  a  of 
this  report  is  the  supplement  of  the  a  of  Ref.  1.  In  addition,  Ca>  Cm, 
and  Cn  are  of  opposite  sign.  (See  the  upper  left  portion  of  Fig.  IV-4 
(Appendix  IV)  for  the  positive  directions  used  in  this  report. ) 
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SECTION  I 
INTRODUCTION 


Initial  aerodynamic  measurements  in  the  reentry  portion  of  the  un¬ 
manned  Apollo  mission  AS -202  trajectory  indicated  several  anomalies 
between  flight  data  and  preflight  predictions  based  on  ground  tests.  In 
an  attempt  to  resolve  some  of  these  anomalies,  a  comprehensive  test 
program  was  conducted  in  the  test  facilities  of  the  von  Karman  Gas 
Dynamics  Facility  (VKF)  (Ref.  1).  As  part  of  this  investigation,  several 
Apollo  models  were  flight  tested  in  the  VKF  1000-ft  Hypervelocity  Range 
(G)  (Armament  Test  Cell,  Hyperballistic  (G)).  The  primary  purpose  of 
these  range  tests  was  to  determine  if  sting  effects  were  significant  in 
the  wind  tunnel  tests. 

Drag  and  stability  measurements  were  made  at  nominal  Mach  num¬ 
bers  of  6.  0  and  8.  5  and  within  a  range  of  free-stream  Reynolds  numbers, 
based  on  the  diameter  of  the  model  (Re^d),  from  50,  000  to  240,  000. 

The  results  of  the  range  tests  are  presented  in  this  report,  and  the 
results  of  the  entire  VKF  investigation  are  summarized  in  Ref.  1. 


SECTION  II 
APPARATUS 


2.1  RANGE 

The  variable-density  range  consists  of  a  10-ft-diam  steel  tank  and 
has  an  840-ft  instrumented  length  (Fig.  1,  Appendix  I).  The  position 
and  attitude  of  the  model  during  the  flight  are  provided  by  43  dual- axis 
spark  shadowgraphs  located  at  nominal  20-ft  intervals  along  the  instru¬ 
mented  range  length.  A  digitized  film  reader  is  used  for  numerical 
interpretation  of  the  shadowgrams.  Angular  orientation  and  linear 
coordinates  of  most  test  configurations  can  be  determined  to  within  ap¬ 
proximately  ±0.  25  deg  and  ±0.  002  ft,  respectively,  at  each  shadow¬ 
graph  station.  A  multichannel  digital  chronograph  records  the  time 
of  model  flight  between  any  two  stations  to  within  ±0.  2  p sec. 

Model  launching  was  accomplished  by  a  two-stage,  powder -hydrogen 
gun  having  a  2.  5 -in.  bore  diameter.  A  blast  chamber  which  absorbs  the 
expanding  muzzle  gases  constitutes  the  initial  85  ft  of  the  range  and  pro¬ 
vides  sufficient  length  along  the  flight  trajectory  for  sabot  separation 
from  the  model  before  entry  into  the  range  proper. 
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2.2  MODELS  AND  SABOTS 

The  configurations  tested  corresponded  to  the  Apollo  command 
module  with  a  smooth  symmetrical  heat  shield  as  shown  in  Fig.  2. 

The  models  were  constructed  of  4340  steel.  Three  different  internal 
model  geometries  were  used  in  order  to  permit  nominal  offset  center- 
of-gravity  (eg)  positions,  e/d,  of  0,  0.  0169,  and  0.  0338  and  correspond¬ 
ing  trim  angles  of  approximately  0,  10,  and  20  deg,  respectively,  in 
flight.  The  fully  offset  eg  position  (e/d  =  0.0338)  was  obtained  with  the 
use  of  two  cylindrical  inserts,  one  of  Kennertium  W2®  and  the  other  of 
LA141A-T7  magnesium-lithium  alloy,  which  were  inserted  into  the 
body  perpendicular  to  the  model  centerline  and  180  deg  apart.  Similarly, 
the  partially  offset  eg  position  was  obtained  by  using  a  single  magnesium- 
lithium  cylindrical  insert.  Each  cylindrical  insert  was  positioned  so  that 
the  extension  of  its  longitudinal  axis  passed  through  the  model  longitudinal 
eg  position;  therefore,  the  three  different  offset  eg  positions  could  be  ob¬ 
tained  without  significantly  affecting  the  longitudinal  eg  position.  Minor 
adjustment  of  the  longitudinal  eg  position  of  a  model  was  made  by  longi¬ 
tudinal  movement  of  a  small  internal  ballast.  Center -of-gravity  loca¬ 
tions  of  the  models  were  measured  to  within  0.  001  in.  in  both  the  trans¬ 
verse  and  longitudinal  directions.  Table  I  of  Appendix  II  presents  a 
summary  of  the  measured  physical  characteristics  of  the  models. 

Four -piece  sabots,  made  of  Lexan®,  were  used  to  position  the 
models  in  the  approximate  trim  attitude  during  the  launching.  Figure  3 
shows  the  angular  alignment  of  a  full -offset  model  in  a  portion  of  its 
sabot.  In  this  photograph  the  model  has  been  rotated  90  deg  about  the 
longitudinal  axis  with  respect  to  its  launch  position  in  the  sabot  in  order 
to  show  more  clearly  one  of  the  inserts  in  the  model. 

In  order  to  determine  roll  rate  during  flight  independent  of  the  con¬ 
ventional  curve -fitting  technique,  described  in  Appendixes  III  and  IV, 
two  small  pins  (not  shown  in  Fig.  3)  were  inserted  into  the  model  apex 
(base)  parallel  to  the  model  centerline  and  approximately  1/8  in.  apart. 
Roll  orientation  as  a  function  of  downrange  distance  traveled  was  then 
obtained  by  determining  the  relative  pin  positions  from  shadowgram 
measurements . 


SECTION  III 
PROCEDURE 


3.1  LAUNCHING 

In  range  testing,  the  transverse  movement  of  a  free -flight  model 
relative  to  the  range  centerline  is  of  concern  in  ensuring  that  the  model 
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remains  within  the  view  area  of  the  shadowgraph  system  during  the 
flight.  The  problem  of  keeping  the  model  within  the  view  area  is  mag¬ 
nified  for  models  having  offset  eg  positions  and  nonzero  lift  coefficients 
at  the  trim  condition.  The  "heavy  model"  type  of  construction  (steel 
models)  used  in  this  investigation  was  a  significant  aid  in  minimizing 
this  transverse  movement. 

The  sabots  used  in  the  present  investigation  were  designed  so  that 
the  models  with  offset  eg  positions  were  oriented  in  the  launcher  at  an 
attitude  approximately  the  same  as  the  anticipated  trim  angle.  Ideally 
only  small -amplitude  angular  motion  would  result,  and  the  trim  angle 
of  a  model  in  free  flight  could  be  determined  from  the  measured  model 
motion  quite  well.  However,  in  the  actual  launchings  of  the  models, 
initial  disturbances  encountered  caused  large -amplitude  motion  in  pitch 
and  yaw  about  the  trim  angle.  Roll  measurements  obtained  with  use  of 
the  pins  previously  discussed  indicated  that  the  roll  rate  of  each  model 
was  approximately  constant  during  its  flight. 


3.2  DATA  REDUCTION 

The  range  data  were  initially  analyzed  using  the  conventional  linear 
analysis  procedures  of  Appendix  III.  No  stability  data  are  presented 
for  models  with  offset  eg  positions  because  it  is  felt  that  violations  of 
equation  restrictions  (discussed  in  Appendix  III)  related  to  the  combina¬ 
tion  of  large  asymmetries  and  large -amplitude  motion  involved  were 
too  severe  to  permit  a  feasible  determination  of  the  stability  derivatives 
for  these  models. 

Typical  angular  motion  plots  are  given  in  Fig.  4  for  representative 
cases  for  each  of  the  eg  positions.  It  is  shown  that  the  assumed  model 
for  the  yaw  program  (Appendix  III)  and  the  values  determined  by  the 
curve  fitting  technique  give  results  that  match  the  experimental  data 
very  closely.  However,  the  values  of  the  trim  angle  determined  by  the 
program  were  found  to  differ  from  the  values  based  on  wind  tunnel  data. 
This  discrepancy  between  the  wind  tunnel  and  range  results  was  found 
to  be  caused  by  the  effect  of  roll  rate  on  trim  angle.  The  magnitude  of 
this  effect  is  given  by  Eq.  (Ill -3)  of  Appendix  III  and  Eq.  (IV- 31)  of 
Appendix  IV.  For  more  typical  configurations  with  higher  fineness 
ratios  the  difference  between  the  rolling  trim  angle,  aq,  an^  the  non¬ 
rolling  trim  angle,  a\Ql  tends  to  be  small,  with  <aq0  smaller  than  a q. 
However,  in  the  case  of  the  Apollo  configuration,  the  inertia  ratio 
(Ix/Iy)  is  greater  than  1  and  aiQ  can  be  larger  than  aq-  Further,  it 
follows  from  Eq.  (Ill- 3)  that  this  difference  is  larger  for  configurations 
with  large  trim  angles. 
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A  more  complete  analysis  of  the  data  was  made  using  the  second- 
order  technique  described  in  Appendix  IV.  This  technique  takes  into 
account  a  number  of  additional  factors  not  included  in  the  yaw  program. 
These  additional  factors  are  the  nonlinearity  of  the  moment  curve,  the 
effect  of  combined  angles  of  attack  and  sideslip  on  the  moment,  and  the 
effects  of  vehicle  asymmetries. 


The  data  analysis  techniques  used  are  summarized  below. 


Parameter  Symmetrical  Models 


Models  with  Offset  eg 


CD 


“t 


^m/tan  ng 


Drag  program 
(Appendix  III) 

Yaw  program 
(Appendix  III) 
(Linear) 

Yaw  program 
(Appendix  III) 
(Linear) 


Drag  program 
(Appendix  III) 

Yaw  program  (Appendix  III) 
and  second-order  method 
(Appendix  IV) 

Not  obtained. 


CL/tan  ae 


Swerve  program  Not  obtained. 

(Appendix  III) 


The  estimated  maximum  probable  errors  in  the  quantities  derived 
from  either  the  linear  or  second-order  analyses  are  listed  below. 


Cd 

±1.5  percent 


j*t 

±1.  25  deg 


Cm/tan  ac 


Cl /tan  a 


±4  percent  ±5  percent 


3.3  DESCRIPTION  OF  SHOTS 

The  physical  characteristics  of  the  models  are  summarized  in 
Table  I,  and  the  test  conditions  and  aerodynamic  data  are  summarized 
in  Table  II.  There  were  three  series  of  shots.  The  first  series  was 
made  with  symmetrical  models  (eg  on  centerline)  at  a  nominal  Mach 
number  of  6.  These  shots  were  numbers  600,  601,  and  602,  Note  that 
the  axial  location  of  the  eg,  x/d,  for  these  shots  was  nominally  0.  242, 
rather  than  0.  270  which  was  the  value  for  the  remainder  of  the  shots. 

The  pitching -moment  coefficient  data  from  shots  601  and  602  given  in 
Table  II  were  adjusted  to  a  eg  location  of  0.  270  for  comparison  with 
other  data.  The  second  series  of  shots  (614  to  620)  was  made  at  a 
nominal  Mach  number  of  8.  5,  at  a  range  pressure  of  from  4  to  8  mm  Hg, 
with  symmetrical  models,  models  with  eg  offsets  of  half  the  design  value, 
and  models  with  full-offset  eg.  The  last  series  included  shots  631,  632, 
and  633;  the  first  two  of  these  were  made  at  a  higher  range  pressure. 
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SECTION  IV 

RESULTS  AND  DISCUSSION 


The  free-flight  range  measurements,  including  trim  angle,  drag 
coefficient,  and  stability  data  are  summarized  in  Table  II.  The  data 
are  presented  in  Figs.  5  through  8,  which  also  include  comparisons 
with  wind  tunnel  data. 

In  Fig.  5a,  range  trim  angle  values  for  the  free-flight  models  in¬ 
cluding  roll  effects,  and  trim  values  corresponding  to  zero  roll,  as 
evaluated  using  the  linear  analysis  of  Appendix  III,  are  presented.  Also 
shown  are  corresponding  wind  tunnel  trim  angle  data.  It  is  apparent 
from  the  range  data  that  roll  effects  on  the  trim  angle  of  an  Apollo -type 
configuration  can  be  significant.  Wind  tunnel  trim  angle  values  (cor¬ 
responding  to  nonrolling  motion)  agree  well  with  the  range  trim  angle 
measurements  adjusted  to  zero  roll  on  the  basis  of  the  conventional 
linear  analysis  of  Appendix  III. 

The  range  trim  angle  data  reduced  using  the  second-order  method 
of  Appendix  IV  are  shown  in  Fig.  5b.  They  indicate  some  improvement 
in  the  spread  of  the  measurements,  but  present  essentially  identical 
results  for  average  trim  angle  of  a  nonrolling  model  as  obtained  by  use 
of  Eq.  ( III  -  3 )  in  the  conventional  analysis. 

With  the  exception  of  shot  632,  the  range  data  are  within  ±1  deg  of 
the  wind  tunnel  data.  This  spread  is  considered  quite  reasonable,  in 
view  of  the  rather  large -amplitude  motions  involved.  No  satisfactory 
explanation  has  been  found  for  the  apparent  deviation  in  the  trim  angle 
of  shot  632,  although  it  should  be  noted  that  a  0.  001-in.  shift  of  eg  posi¬ 
tion  would  place  that  data  point  in  a  reasonably  satisfactory  position. 

The  drag  coefficients  as  a  function  of  effective  angle  of  attack  are 
presented  in  Fig.  6.  Also  shown  are  some  wind  tunnel  data  obtained  at 
comparable  test  conditions.  At  these  relatively  high  Reynolds  numbers 
(Re2d  >  104)  there  is  no  reason  to  expect  a  wall  temperature  effect  on 
the  drag  coefficient  for  this  blunt  shape,  so  the  two  sets  of  data  should 
be  directly  comparable.  At  a  free -stream  Mach  number  of  6,  the  data 
from  the  wind  tunnel  and  the  data  from  the  range  are  in  good  agreement. 
At  a  free -stream  Mach  number  of  8.  5,  the  drag  coefficients  from  the 
range  are  higher  than  those  from  the  wind  tunnel,  especially  for  the 
higher  mean  angles  of  attack.  The  reason  for  this  difference  is  not 
apparent.  The  wind  tunnel  data  indicate  no  Mach  number  effects  on  the 
drag  coefficient  between  free-stream  Mach  numbers  of  6  and  8.  5,  where¬ 
as  the  range  data  indicate  a  slight  increase.  The  differences  between  the 
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two  sets  of  data  are  about  3  percent,  which  is  not  a  large  difference,  but 
it  is  larger  than  the  probable  errors  for  either  set  of  data.  Therefore, 
this  difference  should  be  considered  as  an  indication  of  either  a  slight 
sting  effect  or  a  difficulty  in  obtaining  drag  coefficients  for  large - 
amplitude  motions  at  the  higher  Mach  numbers. 

The  experimental  static  stability  parameter,  Cm/tan  ae,  obtained 
for  the  symmetrical  configurations  using  the  linear  analysis  of  Appendix  III 
is  presented  in  Fig.  7a.  The  values  plotted  correspond  to  measurements 
listed  in  Table  II  that  have  been  adjusted  to  a  longitudinal  reference  point 
at  x/d  =  0.  270.  The  data  indicate  that  the  configuration  is  statically  stable 
at  both  nominal  Mach  numbers  at  which  testing  was  accomplished,  and  no 
Mach  number  effect  is  discernible. 

The  Cp, /tan  ae  values  obtained  for  the  symmetrical  configuration 
are  presented  in  Fig.  7b.  The  data  indicate  negative  values  at  both  Mach 
numbers . 

Linear  curves  with  slopes  corresponding  to  the  reduced  range 
pitching -moment  and  lift-force  derivatives  resulting  from  the  analysis 
of  Appendix  III  are  shown  for  the  M^  =  8.  5  range  data  in  Fig.  8  along 
with  wind  tunnel  Cm  and  Cl  data.  The  wind  tunnel  data  have  been  plotted 
as  a  function  of  tan  a  which  is  more  appropriate  than  a  as  a  parameter  in 
making  the  comparisons  with  the  range  measurements.  Reasonable  agree¬ 
ment  is  indicated  between  the  range  and  wind  tunnel  measurements,  con¬ 
sidering  that  the  range  measurements  are  for  an  amplitude  range  cor¬ 
responding  to  tan  a  up  to  about  0.1.  It  should  be  noted  that  the  decrease 
in  the  absolute  magnitude  of  the  parameter  Cl /tan  a e  with  the  mean  total 
angle  of  attack  as  shown  in  Fig.  7b  is  consistent  with  the  shape  of  the 
wind  tunnel  curve  shown  in  Fig.  8b. 

Some  rough  measurements  of  the  dynamic  stability  damping  param¬ 
eter  were  obtained  for  the  symmetrical  configuration.  The  measurements 
are  not  presented  because  their  questionable  accuracy  does  not  permit 
attaching  more  significance  to  the  measurements  than  making  the  remark 
that  they  indicate  near  neutral  dynamic  stability  for  the  Apollo  configura¬ 
tion. 


It  should  be  noted  that  no  Reynolds  number  effect  is  discernible  in 
the  limited  aerodynamic  data  of  this  report.  The  effect  of  Reynolds 
number  in  wind  tunnel  tests  of  the  Apollo  configuration  is  discussed  in 
detail  in  Ref.  1. 


6 


AEDC-TR- 68-224 


SECTION  V 

CONCLUDING  REMARKS 


The  basic  Apollo  command  module  configuration  has  been  tested  in 
the  VKF  1000-ft  hypervelocity  range  (G)  at  nominal  Mach  numbers  of  6 
and  8.  5  and  within  a  range  of  Reynolds  number  (based  on  model  diam¬ 
eter  and  free -stream  conditions)  from  50,  000  to  240,  000.  Comparison 
of  the  results  with  data  from  wind  tunnels  shows  reasonable  agreement 
with  the  only  difference  worth  noting  being  that  the  drag  coefficient  from 
the  range  tests  is  about  3  percent  higher  than  that  from  the  wind  tunnel 
for  high  amplitude  motions  at  Mm  =  8.5.  Within  this  extent,  sting  effects 
were  found  to  be  negligible. 

The  free-flight  range  experiments  with  an  offset  eg  were  interesting 
in  that  model  roll  in  combination  with  a  large  inertia  ratio  ( (Ix/Iy)  >  1) 
decreased  the  trim  angle  of  the  rolling  Apollo  model  significantly. 
Second-order  effects  on  the  trim  angle  were  suspected  because  of  the 
nonlinear  moment  characteristics  of  the  Apollo  configuration;  however, 
trim  angle  data  on  the  basis  of  either  the  linear  yaw  analysis  or  a  newly 
devised  second-order  analysis  indicate  similar  results. 
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Fig.  6  Drag  Coefficient  as  a  Function  of  Mean  Total  Angle  of  Attack 
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Fig.  8  Comparison  of  Range  and  Wind  Tunnel  Stability  Data  at  M  =  8 
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TABLE  I 

PHYSICAL  CHARACTERISTICS  OF  MODELS 


1.  751 
.  752 


X 

7  d 

0. 

242 

0. 

244 

0. 

242 

0. 

274 

0. 

272 

0. 

273 

0. 

274 

0. 

274 

0. 

274 

0.  270 


0.  270 


0.  0002 


0.  0005 


0.  0004 


0. 0002 


0.  0001 


0.  0157 


0. 0162 


0. 0163 


0. 0164 
0.  0344 


0. 0351 


m  x  102, 
slugs 


1.  282 


1.  279 


1.  361 


1.  357 


1.  279 


1.  284 


1.  286 
1.  364 


1.  360 


Ix  x  10b, 
slug-ft^ 


2.  387 


2.  387 


2.  388 


2.  384 


2.  380 


2.  252 


2.  226 


2.  237 


2.  243 


2.  423 


2.  434 


Iy  x  10b3 
slug -ft  2 


1.  751 


1.  764 


1.  757 


1.  961 


1.  963 


1.  801 


1.  803 


1.  807 


1.  800 


1.  992 


1.  989 


Iz  x  10  , 
slug-ft2 


0.  271 

0. 0339 

1.  358 

2.  407 

1.  985 

0.  275 

0.  0332 

1.  345 

2.  416 

1.  974 

1.  934 


1.  935 


1.  935 


1. 951 


1.  964 


1.  965 


1.  963 
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TABLE  II 

TEST  CONDITIONS  AND  AERODYNAMIC  DATA 


Shot  DRange 

„  ,  Pressure, 

Number  _T 

mm  Hg 


Re„d  x  10  4  Re2d  x  10"4 


P,  6, 

deg /ft  deg 


Drag  Swerve 
Program  Program 


Yaw 

Program 


Second-Order 

Analysis 


Symmetrical 


10.  3 

5.  93 

8.  00 

2.  38 

14.  3 

1.  37 

0.  1 

15.  2 

5-47 

10.  80 

3.  35 

17.  4 

1.  33 

-1.  06 

0.  1 

-0.  126 

31.  3 

5.88 

24.  10 

7.  15 

7.  1 

1.44 

-1.  17 

0.5 

-0.  119 

4.  7 

8.  52 

5.  19 

1.  23 

0.  53 

6.  0 

1.46 

-1.31 

0.  7 

-0.  120 

15.  6 

8.  47 

17.  20 

3.  94 

0.  11 

2.6 

1.  49 

-1.  38 

1.  0 

-0.  118 

Half  Offset 


615 

4.  9 

8.  57 

5.50 

1.  30 

0.  44 

10.  2 

1.  451 

10.4 

11.  5 

10.4 

11.  3 

615 

9.  3 

1.  432 

617 

4.  7 

8.46 

5.  20 

1.  23 

0.  30 

13.4 

1.423 

10.5 

11.0 

10.  6 

11.  0 

617 

13.7 

1.  392 

619 

8.4 

8.  16 

8.  96 

2.  11 

0.  30 

24.1 

10.  1 

10.4 

10.  2 

10.  4 

632 

14.  2 

8.  66 

16.  00 

3.  59 

0.  14 

13.5 

1.  42 

9.7 

9.  7 

9.5 

9.5 

Full  Offset 


4.  8 

8.  37 

5.  24 

1.  25 

0.  34 

15.  5 

1.  39 

21.  2 

23.  3 

21.  0 

4.  7 

8.  72 

5.  37 

1.  24 

0.  47 

11.  8 

1.43 

19.  1 

22.  6 

19.  8 

8.  3 

8.  38 

9.  09 

2.  12 

0.  32 

22.  6 

20.4 

21.  0 

20.  4 

8.  1 

8.  19 

8.  62 

2.  03 

-0.  28 

30.  2 

_ 

19.4 

20.4 

20.  2 

^Determined  Over  the  Interval  from  0  to  600  ft 
^Determined  Over  the  Interval  from  320  to  600  ft 


^Determined  Over  the  Interval  from  0  to  840  ft 
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APPENDIX  III 

BASIC  DATA  REDUCTION  PROCEDURES  FOR  RANGE  TESTS 

This  appendix  describes  the  conventional  techniques  used  in  reducing 
aerodynamic  data  from  measurements  of  the  model  motion  during  the 
test  flight  of  the  model.  Three  different  types  of  data  reduction  programs 
are  involved: 

1.  A  yaw  program  for  determining  the  static  stability  parameter, 
the  damping  parameter,  and  the  trim  angle  from  the  yawing 
motion  of  the  model. 

2.  A  swerve  program  for  determining  the  lift-force  derivative 
from  the  transverse  motion  of  the  model  eg. 

3.  A  drag  program  for  determining  the  drag  coefficient  from  the 
longitudinal  motion  of  the  model. 

1.  YAW  PROGRAM 

In  Ref.  2,  Murphy  derived  a  differential  equation  for  the  motion  of 
an  axisymmetric  model  in  free  flight.  This  equation,  referenced  to  a 
right-hand,  nonrolling  axis  system,  can  be  written  as 

r  +  (h,  -  —  -  ip,) r -  <m  +  iPt -  o  dii-i) 

in  which 

t  V !  .  Wj  .  0 

£  =  -Tj—  +  1  ——  =  sin  Pr  +  1  sin  ar 

v  oo  V  oo 

Additional  quantities  are: 


(See  Fig.  IV- 7  and  the  list  of  symbols  for  definitions  of  the  above 
parameters .  ) 
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It  is  to  be  noted  that  the  lift  and  moment  coefficients  appear  in  the 
terms 

CL  „„j  Cm 

-  and  - 

tan  aE  tan  ae 

rather  than  the  more  conventional  forms 

yCha  and  yCma 

This  has  been  done  to  make  clear  that  the  equation  is  valid  even  for 
large  angles.  It  should  also  be  pointed  out  that  the  coefficients  Cl  and 
Cm  are  to  be  evaluated  at  the  total  angle  of  attack,  a e,  and  need  not  be 
linear  functions  of  angle  of  attack. 

Equation  (III - 1 )  is  valid  for  a  symmetrical  vehicle,  and  in  the  com¬ 
puter  program  the  usual  assumptions  that  y'  =  CLa  =  0  were  made  and 
the  coefficients  were  assumed  constants.  Therefore,  the  program  yields 
effective  mean  values  of  the  quantity  Cm/tan  ae. 

Small  configurational  asymmetries  can  be  handled  approximately  by 
introducing  an  exponential  forcing  function  into  the  right-hand  side  of 
Eq.  (Ill-  1)  (Ref.  2),  and  the  solution  of  the  differential  equation  with 
this  additional  trim  term  can  be  written  as 

£  =  Ai  exp  (fil  +  i  <^/)x  +  A2  exp +  i  02Ox  +  A3  exp  (i  Px)  ( III  —  2 ) 

The  use  of  Eq.  (Ill- 2)  requires  assumptions  of  linear  variation  of  aero¬ 
dynamic  forces  and  moments  with  tan  ae,  equal  static  stability  deriva¬ 
tives  in  pitch  and  yaw,  small  configurational  asymmetries,  and  a 
constant  roll  rate.  Further,  for  small  amplitude  motion  the  assumption 
that  7=1  can  be  made. 

In  a  linear  analysis  of  a  model's  angular  motion,  Eq.  (Ill- 2)  is 
fitted  to  the  measured  components  of  the  complex  yaw  angle,  5,  as 
functions  of  the  distance  traveled,  x,  using  a  differential  corrections 
procedure  (Ref.  3).  There  are  ten  equation  constants  to  be  determined 
in  fitting  Eq.  (Ill -2).  The  A's,  in  general,  are  complex  numbers;  how¬ 
ever,  P  can  be  expressed  as  a  function  of  0  ^  and  0  2-  The  quantities 
and  A 2  depend  on  the  initial  launch  disturbances  encountered  by  the 
model,  but  the  magnitude  of  A3  is  related  to  the  rolling-model  trim 
angle  of  attack  by  the  relation 

sin  at  =  A3 

The  trim  angle  of  a  free -flight  model  corresponding  to  zero  roll  (see 
Ref.  2)  can  be  expressed  as 

ot  =  at  [  1  +  P2  (1  -  Ix/Iy)  ]  (III -  3) 

0  L  J 

Zero  damping,  assumed  in  this  expression,  is  a  quite  valid  assumption. 
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The  stability  parameters  can  be  expressed  in  terms  of  the  deter¬ 
mined  equation  constants  by  the  following  relationships: 


in  which 


tan  ae 


(III-4) 
( HI -  5 ) 


D  =  CD 


— +  kr2  (cm 

tan  ae  \ 


+ 

q 


It  is  apparent  that  if  Eq.  (Ill - 2)  (linear)  is  applied  to  a  model  that 
has  nonlinear  characteristics  the  reduced  stability  parameters  are 
mean  values  corresponding  to  the  motion  interval  involved. 


2.  SWERVE  PROGRAM 

The  lift -curve  slope  is  obtained  from  an  analysis  of  the  transverse 
motion  of  the  model  (Ref.  2).  The  differential  equation  defining  the 
transverse  motion,  referenced  to  a  right-hand,  range-fixed  coordinate 
system  with  the  x-axis  along  the  range  axis,  can  be  written  as 

y'r *  -  -  £  [(drt +  iPdC>0f-  +  (-c?.  -  iC*.> »*p (ms) 

in  which 


Assuming  that  the  coefficients  are  constant,  Eq.  (Ill - 6 )  may  be  inte¬ 
grated  twice  to  yield 

y  +  iz  =  (y0  +  iz0)  +  (yo  +  iz£)x  -  (  Cl  +  iPdCw  ^  (G  +  iH) 

\tan  ae  Pay 

+  (Cyo  +  iCZ())  (I  +  ij)  (HI  -7) 

where 

n  A  *  x2 

G  +  iH  =  — —  r  f  dx,  dx2 

2m  J  J 

0  0 

I  +  ij  =  [(1  -  cos  Px)  +  i(Px  -  sin  Px)] 

2m  P2 

and  an  o  subscript  indicates  the  initial  value.  Equation  (III -7)  is  fitted 
to  the  measured  (y  +  iz)  values  by  a  least-squares  curve-fitting  pro¬ 
cedure  to  obtain  a  value  for  Ci_,/tan  ae. 
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3.  DRAG  PROGRAM 


The  drag  coefficient,  Cp),  was  obtained  from  the  relation 


(III-8) 


in  which  dV/dx  is  the  mean  slope  of  the  curve  of  model  velocity  as  a 
function  of  distance  traveled,  and  Vm  is  the  mean  model  velocity  during 
the  interval  of  the  flight  trajectory  being  examined.  Equation  ( III  -  8 )  is 
based  on  the  assumption  that  Cp)  is  constant  during  the  flight  interval  of 
concern  and  has  been  found  to  provide  quite  adequate  results  when  the 
total  velocity  drop  during  the  flight  interval  is  less  than  about  4  percent 
of  the  initial  velocity.  Since  a  model's  angular  amplitude,  in  general, 
will  vary  along  the  trajectory,  the  drag  coefficient  obtained  by  this  pro¬ 
cedure  corresponds  to  the  mean  total  angle  of  attack  experienced  by  the 
model.  This  mean  angle  of  attack  is  defined  by  the  relation 


-I—  J  [sin  1  x/TiF]  = 


"Li  I 
Li 


(III  -  9) 


where  L2  -  Li  is  the  length  of  the  flight  interval  over  which  Cp>  is  com¬ 
puted. 
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APPENDIX  IV 
SECOND-ORDER  THEORY 


The  comparison  of  trim  angles  based  on  wind  tunnel  data  with  those 
determined  from  the  range  data  by  the  yaw  program  (Appendix  III)  shows 
a  scatter  of  approximately  ±1.  5  deg.  The  second-order  treatment  of 
this  appendix  was  undertaken  in  an  attempt  to  reduce  this  scatter  by 
accounting  for  certain  known  small  effects  with  the  aid  of  available  wind 
tunnel  data.  The  result  leads  to  only  a  minor  reduction  in  the  scatter. 
Nevertheless,  the  complete  treatment  is  included  because  it  points  out 
some  interesting  peculiarities  in  the  moment  curve  of  the  particular 
shape  being  studied,  and  it  presents  in  a  specific  example  an  extension 
of  the  basic  work  of  Murphy  (Ref.  2)  to  the  case  of  a  vehicle  with  asym¬ 
metries  which  are  not  necessarily  small. 


DYNAMIC  EQUATIONS 


The  yaw  program  (Appendix  III)  is  written  in  terms  of  the  range 
axis  system,  using  the  angles  ar  and  /3r.  This  is  the  most  advantageous 
system  for  a  completely  symmetrical  vehicle.  During  the  tests  covered 
in  this  report,  conditions  were  encountered  for  which  the  assumption  of 
symmetry  was  not  justified,  so  that  it  was  necessary  to  obtain  results 
with  respect  to  a  body-fixed  axis  system. 


The  following  equations  for  a  body- fixed  axis  system  with  products 
of  inertia  equal  to  zero  are  from  Etkin  (Ref.  4). 


M  =  !y  -  FIy  PR  =  Cm  qAd 

dt  1 

N  =  Iz  -iB.  +  GIZ  PQ  =  Cn  qAd 

dt 

in  which 


(IV- 1) 
(IV-2) 

(IV-3) 

(IV-4) 


The  body-fixed  axis  system  is  shown  in  Fig.  IV- 1.  The  angular  veloci¬ 
ties  are  in  units  of  radians  per  second. 
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The  relationships  among  the  velocity  components  and  the  angular 
velocity  components  are  also  obtained  from  Etkin. 

y—  +  RU-PW  =  Y/m  +  g  cos  Q  sin  $  (IV-5) 

—  +  PV-QU  =  Z/m  +  geos  6  cos  $  (IV-6) 

dt 

in  which  Y  and  Z  are  the  aerodynamic  forces,  and  6  and  $  are  the  eleva¬ 
tion  and  roll  angle,  respectively,  of  the  axis  system  relative  to  an 
inertial  system.  The  approximation  is  now  made  that  the  right-hand 
sides  of  these  two  expressions  may  be  neglected.  Neglecting  the  gravi¬ 
tational  terms  can  be  easily  justified.  The  aerodynamic  force  terms, 
if  included,  would  appear  as  damping  terms  in  the  final  equations.  Pro¬ 
vision  for  damping  will  be  added  at  a  later  stage  in  the  development. 
Neglecting  the  damping  terms  at  this  point  is  equivalent  to  neglecting  the 
differences  in  the  damped  and  undamped  natural  frequencies,  which  can 
be  justified  on  the  basis  of  the  small  damping  associated  with  the  partic¬ 
ular  configuration  being  considered. 


The  velocity  components  V  and  W  are  now  eliminated  by  using  the 
angles  of  attack,  a,  and  sideslip,  /3,  defined  by  (Fig.  IV-2) 


U  =  Voo  cos  a  cos  /3 
V  =  Voo  cos  a  sin  /3 
W  =  Voo  sin  a 

Voo  =  yj  U2  +  V2  +  r 


(IV-7) 


Including  only  terms  up  to  second -order  in  a  and  j3  simplifies  these  to 

U  =  Voc  [1  -  («2  +  £2)/2] 

V  =  VacjS 

W  =  Vooa 


These  results  may  now  be  used  in  Eqs.  (IV-5)  and  (IV-6).  In  doing  so, 
P  is  considered  a  quantity  of  zero  order,  and  Q  and  R  are  considered 
quantities  of  first  order.  At  the  same  time,  the  conversion  is  made 
from  time,  t,  to  distance,  x,  as  the  independent  variable,  by  using 


dx 

dt 


v, 


The  angular  velocities  are  then  measured  in  radians  per  foot.  The 
results,  correct  to  the  second  order,  are 


+  r  _  Pa  =  o  (IV-8) 

-jjs  +  P£  -  Q  =  0  (IV-9) 
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These  are  now  used  in  Eqs.  (IV- 1)  and  (IV- 2)  to  give 
— - FP2a  +  P(1  +  F)  —  =  Cm  qAd- 

dx2  dx  V2„o  Iy 


(IV-10) 


-P(l  +  G) 


-  GP2/3  =  - 


cn  qAd 
V2~  I, 


(IV-11) 


Further  derivation  depends  upon  determining  the  functional  dependence 
of  the  moment  coefficients  on  the  angles  a  and  /3. 


MOMENT  COEFFICIENTS 

Figure  IV- 3  gives  wind  tunnel  data  for  the  pitching-moment  and 
axial -force  coefficients  as  a  function  of  angle  of  attack  from  Ref.  1. 
These  may  be  represented  by* 

Cmi  =  -0. 00220a  +  0.0000085  a2  (IV-12) 

Ca  =  1.465  -  0.000339  a2  (IV- 13) 


The  pitching  moment  Cmi  is  taken  about  the  point  on  the  axis  of  the 
vehicle  at  which  the  plane  perpendicular  to  the  axis  which  contains  the 
offset  eg  intersects  the  axis  (Fig.  IV-4).  The  pitching  moment  about 
the  offset  eg  is 


Let 


Cm  =  CA  |  +  Cmi  (IV- 14) 

=  0.03384  P0 
d 


in  which  PQ  represents  the  fraction  of  offset  (i.  e.,  PQ  =  1  corresponds 
to  the  full-offset  case).  Then 


Cm  =  (1.465  -  0.000339  a2)  (0.03384  P0)  -  0.00220  a  +  0.0000085  a2  (IV-15) 


*The  data  given  in  Fig.  IV-3  are  more  accurately  represented  by 

cmi  =  -0.00208a  +  0.0000050  a2 

for  the  smaller  angles.  The  representation  given  by  Eq.  (IV-12)  was 
used  in  the  second-order  data  reduction  program  because  it  represents 
the  data  better  at  the  large  angles.  The  effects  of  the  slight  differences 
are  negligible  for  the  development  in  this  appendix. 
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From  this  relation,  the  trim  condition  can  be  determined.  The  follow¬ 
ing  results  are  obtained: 

PQ  <*t>  de§ 

0  0 

0.5  11.4 

1.0  21.9 


The  representation  for  the  pitching  moment  given  by  Eq.  (IV- 15)  is 
valid  only  if  the  velocity  vector  lies  in  the  x-y  plane.  If  the  velocity 
vector  is  outside  the  x-y  plane  (i.  e. ,  the  vehicle  is  at  a  sideslip),  then 
more  care  must  be  taken  in  choosing  the  representation. 


The  location  of  the  eg  is  unimportant  for  determining  the  magnitude 
of  the  aerodynamic  moment  and  C^.  Therefore,  the  magnitudes  of 
these  quantities  are  determined  by  the  total  angle  of  attack,  ae,  which 
is  defined  as  the  angle  between  the  x-axis,  or  the  axis  of  symmetry, 
and  the  velocity  vector,  measured  in  the  plane  containing  the  two. 
Therefore,  if 


=  \}~a2  +  P2 


CA  =  1.465  -  0.000339 


and 


'll 


=  -  0.00220  ae  +  0.0000085  ae2 


(IV-  16) 
(IV- 17) 

(IV- 18) 


The  projections  of  the  total  moment  coefficient  on  the  appropriate  axis 
thqn  give  the  pitching-  and  yawing-moment  coefficients  relative  to  the 
axis  system  fixed  in  the  offset  eg  (Fig.  IV-4) 


Therefore, 


in  which 


a/ ae 

Cn  =  Cm ^ ^  /3/ ae 

Cm  =  -0.00220  a  +  C2  aae  +  0.0496  P0  -  ae 
Cn  =  0.00220  /3  -  C2/3ae 

C,  =  0.0000115  P0 
C2  =  0.0000085 


(IV  - 19) 
(IV-20) 

(IV-21) 

(IV-22) 

(IV-23) 

(IV-24) 


Figure  IV-5  shows  the  resulting  pitching  moment  as  a  function  of  a  and  )3, 
and  shows  that  a  vehicle  undergoing  a  yawing  oscillation  will  trim  at  a 
lower  value  of  a  than  one  that  is  not  yawing. 
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FURTHER  DEVELOPMENT  OF  EQUATIONS 


The  expressions  for  the  pitching  and  yawing  moments  are  now 
introduced  into  Eqs.  (IV- 10)  and  (IV- 11).  The  damping  terms  will  be 
added  later. 

+  ( 0.00220  ■  -q-A-~  -  FP2)  a  +  P(1  +  F) 

dK  V  v~  Iy  / 

=  qAd--  (C2aae  +  0.0496  P0  -  C,  ae2) 

Vie  Iy 

-P(l  +  G)  —  +  -2-2-  +  (0.00220  — Shi  -  GP2)  /3  = 

dx  dx2  \  V^o  Iz  / 

The  following  development  has  two  objectives: 

1.  The  equations  are  to  be  written  with  all  linear  terms  on  the 
left  side  and  the  nonlinear  terms  on  the  right-hand  side. 

2.  The  numbers  0.  00220  and  0.  0496  PD  appearing  in  the  moment 
expressions  will  be  eliminated  in  favor  of  the  natural  frequen¬ 
cies  and  the  trim  angle,  which  are  the  quantities  that  may  be 
determined  from  the  range  tests. 

The  nonlinear  quantities  and  C2  will  be  retained  at  the  values  given 
by  the  wind  tunnel  tests,  since  it  is  thought  that  it  is  not  possible  to 
evaluate  these  small  parameters  from  the  range  data. 


_d£_ 

dx 


(IV-25) 


-241.  C2jS «e  (IV- 26) 

v  2  T 

V OO  Az 


From  Eqs.  (IV-25)  and  (IV- 26),  the  steady-state  condition  is 


/3  =  0 

a.  =  Q-e  =  o.  t 


(IV-27) 


so 


0.00220 


qAd 
Voo2  Iy 


at 


-3^i-  [(C,  -  C,)  at2  +  0.0496  P0] 

Voi  Iy 


(IV-28) 


At  this  value  of  a,  we  find  from  Eq.  (IV- 21) 

Cmao  =  -0.00220  +  2(C2  -  C,)  at  (IV-29) 

It  is  therefore  convenient  to  define  the  oscillating  frequency  for  small 
oscillations  by 

cmao  qAtl 

=  "  v J  iy  (IV- 30) 

=  [0.00220  -  2(C2  -  CJ.atl  r£tr~ 
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The  effect  of  the  roll  rate  on  the  static  stiffness  is  apparent  from 
Eq.  (IV-28).  In  order  to  compare  with  wind  tunnel  data,  the  trim  angle 
for  zero  roll  rate  must  be  determined,  whereas  the  range  data  "solve" 
Eq.  (IV-28)  with  non-zero  roll  rate.  Since  the  correction  for  roll  rate 
is  small,  the  second-order  terms  of  Eq.  (IV-28)  (those  involving  C2  ~C\) 
may  be  neglected  in  evaluating  the  correction.  With  this  assumption,  the 
trim  angle  without  roll,  aiQ,  is  related  to  that  with  roll,  a t,  by 

ato  =  at[l  -  F  (-£)]  (IV- 31) 

The  parameters  P,  and  ua  are  to  be  determined  from  the  range  data 
by  the  technique  described  later. 


In  order  to  eliminate  the  parameters  obtained  from  the  wind  tunnel 
data,  the  term  0.  0496  PQ  in  Eq.  (IV-25)  is  eliminated  by  using  Eq.  (IV-28). 
The  term  involving  the  parameter  0.  00220  is  eliminated  by  using 
Eq.  (IV- 30).  In  addition,  the  following  substitutions  are  made: 

a  =  at  +  a,  (IV-32) 


K,  =  Ct/\Cma  | 

K2  =  C2/|Cm  | 
w’0 

Equation  (IV-25)  then  becomes 

+  fe,  *  _  FP2)  a,  +  P(1  +  F)  ^  =  coa2  [atae  +  -  at2 

dx2  dx 

-  2atai)  -  Kj  (Kl2  +  (S2)] 

Similarly,  for  the  yaw  equation,  let 

=  0.00220  —  C2at 

Cnj3o  qAd 

CD  a  =  - - - 

P  V  2  T 

v  OO  Lz 

K3  =  C2/Cnpo 

Equation  (IV- 26)  then  becomes 


(IV-33) 

(IV-34) 


(IV-  35) 

(IV-  36) 
(IV-37) 
(IV-38) 


-P(l  +  G) 


daj 

dx2 


+  -4  +  (co2a  -  GP2)  |8  = 

dx2  P 


p 


K3  (ae  -  at) 


(IV-39) 


One  further  simplification  is  introduced.  In  obtaining  the  second- 
order  solution,  the  term  ae  is  troublesome.  Therefore,  an  approxima¬ 
tion  will  be  made  to  the  value  of  this  term.  In  general, 

ae  =  yjo-t  +  2a  [t*!  +  a 2  +  /32  (IV-40) 
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Assume  that  is  the  dominant  term,  and  the  remaining  terms  are 
small.  Then 

a\  +  p2 

ae  ~  at  +  (IV-41) 


This  is  used  in  Eqs.  (IV-35)  and  (IV-39)  to  give 


and 


in  which 


d2  a , 
dx2 


+  c2  ax  +  P(1  +  F)  —  =  a>a  OW  +  F2/32] 

A  dx 


-  P  (1  +  G) 


daL 
d  x 


d*<3 

dx2 


cu2b:/8  =  co2^  K 3a,^ 


Ft  =  1.5  Ka  -  K, 
F2  =  0.5  K,  -  K, 


c2a  -  FP2 
CU^g  -  GP2 


(IV-42) 

(IV-43) 


(IV-44) 

(IV-45) 


In  obtaining  these  results,  third-order  terms  in  a\  and  /3  have  been 
neglected. 


LINEAR  SOLUTIONS 


The  linear  solution  to  Eqs.  (IV-42)  and  (IV-43)  is  obtained  by 
neglecting  the  right-hand  sides.  The  resulting  equations  to  be  solved 
are 

d  CL  1  (i  ^ 

— +  <u\  ai  +  P  (1  +  F)  —  =  0  (IV-46) 

~  P(1  +  G)  +  -0-  +  «b  i3  =  0  (IV -4 7) 


The  characteristic  equation  is  given  by 

s2  +  co\  PCI  +  F)s 

A  =0 
-P(l  +  G)s  s2  +  co2B 

or 

s<  +  [<u2b  +  <o2B  +  (1  +  F)(l  +  G)P2]  s2  +  cu2a  «2b  =  0 
This  may  be  factored  to 

(s2  +  a>2)  (s2  +  oi22)  =  0 


(IV-48) 
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The  relations  among 
&>!  =  0.5 
c o2  =  0.5 

co  A  =  0.5 
co  b  =  0.5 


<*>1,  W2> 

r 

^A> 

and 

WB 

are  found  to  be 

1 

+ 

cob )2 

+  z. 

+  S 

r 

■\j  (wA  - 

cob)2  +  ZjJ 

■\J  (a  A 

+ 

«b)2 

+  z, 

-  s 

\j((0  A  - 

cob)2  +  Z, 

^((Ol 

+ 

co2)2 

-  Z, 

+  s 

-\j  (cox  - 

co2)2  -  Z, 

" 

\l  ^ 

+ 

o2y 

-  Zx 

-  s 

V  (w‘  “ 

tsT 

1 

o« 

3 

in  which 


Z,  =  (1  +  F)  (1  +  G)P2 

S  =  +  1  co^  y  co  p  |j  co  i  y  co2 

-  1  CO  X  <  COp‘,  &>!  <  co2 


(IV-49) 


(IV-50) 

(IV-51) 


The  solution  to  the  equations  may  be  written 

=  Da  sin  (&jlx  +  8a)  +  cos  (co2x  +  8p)  (IV-52) 

jS  =  sin  ( co2x  +  8p)  ~  cos  (c^x  +  8a )  ( IV  —  5 3) 

The  solution  has  been  set  up  so  that  as  the  roll  rate  P  goes  to  zero,  the 
solution  with  frequency  becomes  the  oscillation  in  a.  The  parameters 
Ti  and  T2  are  given  by 


Co\  ~  Co\  COfjCOi  —  C0/^C02 

CO t  P(1  +  F)  ~  coB  PG  +  p) 


(IV-54) 


a2  -  m  B  C0^C02  -  COBCOj 

co2  P(1  +  G)  «Ap  (1  +  G) 


1  +  F  fc»B  rp 

1  +  G  co  \ 


(iv-55) 


It  should  be  noted  that  for  a  completely  symmetrical  vehicle,  the 
above  results  yield 

wi  -  co2  =  S(1  +  F)P 
T,  =  S 
t2  =  -S 


PARTIAL  SECOND-ORDER  SOLUTION 

One  method  of  obtaining  a  solution  to  Eqs.  (IV-42)  and  (IV-43)  is  to 
use  the  perturbation  method  (Ref.  5).  The  technique  is  based  upon  the 
assumption  that  the  nonlinear  terms  are  small,  so  that  an  approximate 
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evaluation  of  these  terms  can  be  used.  For  the  particular  problem 
under  consideration,  this  involves  using  the  solutions  to  the  linear  case 
as  given  by  Eqs.  (IV-52)  and  (IV-53)  in  the  right-hand  sides  of 
Eqs.  (IV-42)  and  (IV-43),  and  obtaining  new  solutions  for  a\  and  j3.  The 
solution  will  then  appear  as  the  linear  solution  plus  additional  second- 
order  terms. 


Upon  combining  Eqs.  (IV-42)  and  (IV-43)  for  ai  and  j3,  there  is  ob¬ 
tained 


<P  at 

dx2 


+  W  +  0>22)  +  =  "a2  (-77  +  “  b)  (F.  «.2  +  F2i32)| 

dx  \dx  ]  ' 


-  co2p  K3  P(1  +  F)  («iyS)| 


(IV-56) 


d4/3 


+  (ft)!2  +  ft)22)  +  ft)!2ft)22/3  =  ft)2^  K3  (-E-  +  ft)2A)  («!^)| 

dx  dx  1 


+  P(1  +  G)  ft)«2  -f  (F idi  +  F2/32)| 

ax  n 


(IV-57) 


with  the  understanding  that  the  expressions  for  and  |3  given  by 
Eqs.  (IV-52)  and  (IV-53)  are  used  in  the  right-hand  sides  of  these 
equations. 


It  should  be  obvious  that  the  resulting  solution  is  quite  lengthy.  It 
was  decided  to  limit  the  solution  to  the  two  largest  terms,  the  constant 
term,  and  the  terms  involving  sines  and  cosines  with  frequency  -  U2- 
There  are  also  terms  involving  sines  and  cosines  with  frequencies 
2 ui,  2u>2,  and  +  to 2 •  The  coefficients  of  these  terms  are  small  be¬ 
cause  of  the  higher  frequencies,  and  furthermore,  the  data  from  the 
range  seldom  covered  more  than  one  cycle  of  frequency  co^  or  u>2>  so 
that  these  terms  have  not  been  included  in  the  remaining  development. 


With  these  simplifications,  the  solution  becomes 

Oo  +  “f-  [(Ft  +  T32F2)  Da2  +  (F*  +  T22  FJ  D2  ]  +  Dasin(ft)!x  +  8a) 

2ft)  A  P 

+  T2Dp  cos  (ft)2x  +  )  +  DaDp  [ft)a2  (F3T!  +  F2T3)  (w^  -  (ft)3  -  co2)2) 

,  .  sin  (ft),x  +  Sn  —  C0,x  —  8a 

+  0.5  K3P(1  +  F)  co2g  (ft)t  -  ft)2)(l  -  T/T*)]  - - - - - - - ^ -  /T,r 

p  ft)ift)2  (2ft)t  -  ft)2)  (2ft)2  -  ft)t)  (IV-58) 


fi  —  Dp  sin  (cj2x  +  5p)  —  T^Do;  cos  (ft)tx  +  8a)  +  DaDp  [&)p  K3  (1  —  TiT2) 

(ft)2A  -  (ft)!  -  ft)2)2)/2  +  ft)a2(ft)!  -  ft)2)P(l  +  G)(F!T2  +  F2T>)] 

COS  (ft)! X  +  8a  -  ft)2X  -  Sp  ) 

ft)l&)2  (2 ft)!  -  ft)2)  (2ft)j  -  ft) ! )  (IV-59) 

These  are  the  equations  that  were  fitted  to  the  data. 
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RELATIONSHIPS  INVOLVING  RANGE  ANGLES  AND  BODY-FIXED  ANGLES 

The  data  from  the  range  consist  of  values  of  the  elevation  angle, 
ar,  and  the  yaw  angle,  j3r,  at  various  stations  along  the  range.  Rela¬ 
tionships  must  be  found  between  these  angles  and  angles  a  and  ft  appear¬ 
ing  in  the  equations  derived  in  this  appendix.  Referring  to  Fig.  IV-6, 
two  axes  systems  are  shown,  an  xyz  system,  fixed  in  the  model  (as 
shown  in  Fig.  IV- 1),  and  an  xryrzr  system  fixed  in  the  range,  with  xr 
directed  along  the  range  axis,  zr  vertically  downward,  and  yr  hori¬ 
zontal  and  forming  a  right-hand  system.  The  vector  a  is  the  vector 
from  the  center  of  mass  perpendicular  to  the  axis  of  the  vehicle 
(Fig.  IV- 1).  The  orientation  of  the  xyz  system  relative  to  the  xryrzr 
system  is  specified  by  three  angles  defined  as  follows: 

1.  Beginning  with  the  xryrzr  system,  rotate  about  the  xr  axis 
through  the  roll  angle,  $,  to  the  Xj-y-jz^  axis  system,  which 
contains  the  center  of  mass  and  the  vector  a. 

2.  Rotate  about  the  y^  axis  through  the  angle  of  attack,  a,  to  the 
x^y^z  axis,  so  that  the  vector  a  lies  along  the  z  axis. 

3.  Rotate  about  the  z  axis  through  the  sideslip  angle,  jS,  to  the 
xyz  axis  system,  so  that  the  x  axis  is  parallel  to  the  axis  of 
the  vehicle. 

The  unit  vector  i  along  the  x  axis  of  the  vehicle  can  be  related  to  the  unit 

-4  ->  -» 

vectors  ir,  jr,  and  kr  of  the  xryrzr  system.  The  result  is 

-»  -» 
i  =  cos  a  cos  /3  ir  +  (sin  a  cos  /3  sin  4>  -  sin  /3  cos  <j>)  jr 

-  (sin  a  cos  /3  cos  <jS  +  sin  /3  sin  0)  kr  (IV-60) 

The  range  angles  are  the  elevation  angle,  av>  and  the  yaw  angle,  f3r. 
These  are  defined  with  respect  to  an  x-^y^zj  coordinate  system,  which 
has  its  xi  axis  along  the  x  axis  of  the  vehicle,  and  the  y^  axis  horizontal. 
This  is  therefore  a  nonrolling  body  axes  system.  The  orientation  of  this 
axis  system  relative  to  the  range  axis  system  is  specified  by  two  angles, 
ft}  and  0]_,  defined  as  follows  (Fig.  IV-7). 

1.  Rotate  about  the  zr  axis  through  an  angle  /3-j_  to  the  xjy^zr  axis 
system,  so  that  the  x  axis  lies  in  the  xjzr  plane. 

2.  Rotate  about  the  y^  axis  through  an  angle  6,  to  the  x^y^z-^  axis, 
so  that  the  x^  axis  coincides  with  the  x  axis. 

For  this  transformation 

^  ^  ^  t  i 

i  =  cos  /3i  cos  0,  ir  -  sin  /3i  cos  0,jr  -  sin  6l  kr  (IV- 61) 
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Equating  the  corresponding  coefficients  of  Eqs.  (IV-60)  and  (IV-61) 
gives 

cos  cos  dl  =  cos  a  cos  /3 

sin  /3,  cos  0,  =  -  sin  a  cos  /3  sin  +  sin  /3  cos  (f> 

sin  dl  =  sin  a  cos  /3  cos  <^>  +  sin  /3  sin  <f) 

Therefore, 

tan  /3,  =  -  tan  a  sin  <^>  +  tan  /3  cos  0/cos  a  (IV~6  3) 

tan  0,  =  cos  /3,  [tan  a  cos  0  +  tan  /3  sin  0/cos  a]  (IV-64) 

The  angles  ar  and  /3r  are  defined  such  that  (Fig.  IV -7) 


( 

.  (IV-62) 


V  =  Voo  sin  /Sr 

W  =  Voo  sin  aT 

Therefore, 

sin  =  sin  /3(  I 

sin  ar  =  cos  /3,  sin 


For  small  angles  in  0^  and  /3j,  these  simplify  to 


«r 

/^r 


cos  <j>  a  +  sin 
-  sin  (f>  a  +  cos 


P  l 

<$>  Pj 


(IV-65) 


(IV-66) 


(IV-67) 


Both  the  nonlinear  relations  (Eqs.  (IV-63),  (IV -64),  and  (IV-66))  and  the 
linear  relations  (Eq.  (IV-67))  have  been  used  in  the  data  reduction  pro¬ 
gram,  and  differences  were  noted  for  only  one  run. 


The  roll  angle  <t>  can  be  expanded  by 

<f>  =  Px  +  4>0 


(IV -6  8) 


SOLUTION  PROCEDURE 

The  solutions  given  by  Eqs.  (IV-58)  and  (IV-59)  give  the  angle  of 
attack,  a j,  and  the  sideslip  angle,  ft,  in  terms  of  the  distance  x  and  the 
parameters  ua,  up,  P,  aQ,  D#,  D^,  8a,  and  6p.  In  addition  a  damping 
term  of  the  form  e--^-  was  applied  to  the  first-order  oscillatory  terms. 
The  same  value  of  D  was  used  for  both  modes,  since  the  damping  was 
quite  small. 
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The  solution  procedure  consisted  of  using  the  method  of  differen¬ 
tial  corrections  (Ref.  3)  to  determine  the  parameters  %,  U R,  P,  a Q, 
Dff)  Dj3,  6a,  6 ft,  (j>Q,  and  D  to  minimize  the  sum  of  the  squares  of  the 
differences  between  the  calculated  values  of  ar  and  fir  and  the  experi¬ 
mental  values. 


39 


AEDC-TR-68-224 


Fig.  I V- 1  Body-Fixed  Axi  s  System 
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Fig.  IV-3  Wind  Tunnel  Data  for  Moment  and  Axial  Force 
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Fig.  IV-4  Projection  of  Moments 
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Angle  of  Attack,  a,  deg 

Fig.  IV-S  Pitching-Moment  Coefficient  about  Full-Offset  Center  of  Gravity 
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